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We study screening correlators of quark-antiquark composites at T = 2T C , where T c is the QCD 
phase transition temperature, using overlap quarks in the quenched approximation of lattice QCD. 
As the lattice spacing is changed from 1/4T to a = 1/6T and 1/8T, we find that screening correlators 
change little, in contrast with the situation for other types of lattice fermions. All correlators are 
close to the ideal gas prediction at small separations. The long distance falloff is clearly exponential, 
showing that a parametrization by a single screening length is possible at distances z > 1/T. The 
correlator corresponding to the thermal vector is close to the ideal gas value at all distances, whereas 
that for the thermal scalar deviates at large distances. This is examined through the screening 
lengths and momentum space correlators. There is strong evidence that the screening transfer 
matrix does not have reflection positivity. 
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I. INTRODUCTION 

Current experiments at the Brookhaven RHIC and an- 
other to be started soon at the CERN LHC are engaged 
in creating and understanding the phase of matter called 
the quark gluon plasma. The existence of this phase was 
predicted by lattice QCD, and this technique has been 
used to study many properties of the plasma. In spite of 
this, much remains to be done. One of the major ques- 
tions, that which concerns us here, is the nature of the 
interactions between quasiparticle excitations. 

Decisive information on flavoured quasiparticles has 
come from a recent study of the linkage between quantum 
numbers [l| . A rapid and clear change in the linkage be- 
tween quantum numbers was seen across the QCD phase 
transition. In the high temperature phase flavour quan- 
tum numbers such as baryon number, charge, strangeness 
etc., are linked to each other in ways that suggest that 
the light flavoured excitations at finite temperature QCD 
include quarks. 

At the same time it is clear that the interactions be- 
tween these quasiparticles cannot be small. At length 
scales of order the QCD coupling, g, is of order 1 
0, leading to well known problems. For example, the 
Debye screening length of gluons is more complicated 
than one expects in perturbation theory Q and contains 
pieces which are entirely non-perturbative The non- 
perturbative pieces have been isolated and studied on the 
lattice [5|. Interestingly, when one studies the screen- 
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ing of glueball-likc quantities, i. e., spatial correlations of 
colour singlet operators made out of gauge fields, then a 
surprising simplicity arises above about 1.25T C . Taking 
two quantum number channels which can be constructed 
by a minimum of two and three gluon operators, the ratio 
of the screening masses is seen to be close to 3/2 @. 

There are parallels in the study of hadron-like screen- 
ing correlators in the QCD plasma. Parity doubling 
was seen in the high-temperature phase of QCD in the 
first study of these correlators using staggered quarks [|[ . 
Further, this first study already showed that the screen- 
ing mass from the baryon-like correlator was 3/2 times 
the screening mass from some of the meson-like correla- 
tors. A finite size scaling study using staggered quarks 
Q showed that a particular combination of the vector 
and axial-vector meson-like correlators was very close to 
the free-field theory (i. e., ideal gas) predictions. This 
behaviour is generic, bein g se en in quenched [l(| and 
dynamical QCD with two |ll| and four [l§] flavours of 
staggered quarks as well as with Wilson quarks [l3| . Ex- 
trapolation to the continuum using staggered quarks (l7| 
showed that the remaining screening masses approached 
their ideal-gas values. However, the correlators differed 
strongly from the ideal-gas correlators at small distances. 
This was also seen in a later study of screening correlators 
with Wilson quarks [l8[ . 

Overlap quarks [l|| have the advantage of preserving 
all chiral symmetries on the lattice for any number of 
massless flavors of quarks [2(| • This is in contrast to 
other formulations, such as Wilson's, which break all 
chiral symmetries, or the staggered, which break them 
partially. Since the number of pions and their nature 
is intimately related to the actually realized chiral sym- 
metry on the lattice, one expects any realization of chiral 
quarks on the lattice to provide insight into this question. 
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TABLE I: The break up of irreps under the successive break- 
ings of the symmetries of the transfer matrix. The Af com- 
ponents of the V/AV correspond to the t polarization and the 
E to the x and y polarizations. All states also carry a la- 
bel for charge conjugation, C. This has been dropped in this 
paper since we consider only the C = 1 states. 



In our earlier study of screening correlators using over- 
lap quarks j25| with lattice spacing a — 1/4T, the cor- 
relation functions were found to decay exponentially at 
large separation, and showed none of the fine structure 
that plague staggered and Wilson quarks. As a result, 
the screening mass was a good parametrization of the 
screening correlator at large distances. The screening 
masses in all channels were closer to the expected weak- 
coupling continuum limit for overlap quarks than stag- 
gered and Wilson. However, that from the (would be 
zero-temperature) scalar and pseudoscalar was lower by 
about 10% than the others. 

In this paper we extract complete information on 
meson-like screening correlators by extending the anal- 
ysis of [1, [lU to overlap quarks. At the same time, we 
extend the analysis towards the continuum limit by using 
lattice spacings of a = 1/6T and 1/8T. With these two 
inputs we are able to resolve all the currently outstanding 
questions on screening correlators. Parts of these results 
were presented in [la ]. 



II. SYMMETRIES 

We study correlation functions of colour singlet oper- 
ators constructed with a quark and an antiquark — 



C T (z) 
M T {z) 



(Mr(O)Mt(z)), 



where 
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ip is a Dirac spinor, T is a Dirac matrix, the angular 
brackets denote averaging over gauge configurations and 
the quark bilinear Mr is projected to zero momentum in 
the slice orthogonal to z by the summation. By choosing 
different Dirac matrices, T, one explores different quan- 
tum numbers, and the usual nomenclature is explained 
in Table HI Note that the correlators projected on zero 



momentum {k x = k y = k t = 0) have mass dimension 3, 
and hence the quantity a 3 C(z) is dimensionless. 

Since we work at finite temperature, the Dirac opera- 
tor is defined with antiperiodic boundary conditions on 
t (1 < t < N t ). We chose to impose periodic boundary 
conditions in the spatial directions (1 < x, y < N s and 
1 < z < N z with N s < N z ). Later we shall have occasion 
to use the aspect ratio ( = N z /N s . The sum over x, y 
and t in eq. ((TJ) accomplishes a projection on to zero total 
momentum in these three directions. 

Due to the inequivalence of the spatial directions x and 
y with the Euclidean time direction t at finite tempera- 
ture, a slice of the lattice orthogonal to the z direction 
differ at zero and finite temperature. These are part of 
the symmetries of the transfer matrix. 1 Therefore the 
classification of operators Mp differ at zero and finite 
temperature. This is described below and summarized 
in Table [J 

At T = in the continuum one has the full rotational 
symmetry 0(3) and the discrete symmetries of parity 
P and charge conjugation C. States corresponding to 
the irreducible representations (irreps) of 0(3) x Z 2 (P) x 
Z 2 (C) are labeled by J PC . On the lattice this group 
is broken to the discrete subgroup which is the cubic 
group O h . The PC states (S/PS) become the A PC of 
the cubic group and the 1 states (V/AV) become the 
three dimensional irrep F PC (sometimes called the T PC ) 
[22l ]. The 2 PC breaks into a three dimensional irrep called 
the F PC (or T PC ) and a two dimensional irrep called the 
E PC . 

For T 7^ 0, the continuum theory has a 0(2) rotational 
symmetry in the xy plane and a Z2{T) symmetry (Eu- 
clidean time reversal) in the Euclidean time direction. 
This constitutes the cylinder group, 0{2) x ^(T), which 
is a subgroup of 0(3) x Z 2 (P). On the lattice this is 
broken to the dihedral group D\ = L> 4 x Z 2 (T) (where 
D4 is the group of symmetries of a square), which is also 
a subgroup of Oh, as is to be expected [21j . 

The J irreps of T = break up into irreps 
at finite temperature, where M — J z (since the cylinder 
group is Abelian, its complex irreps are one-dimensional) . 
The reduction under the lattice symmetries is shown in 
Table [I] The group theory is the most general. How- 
ever, it does not restrict special cases in which some of 
the independent representations may become identical 
for dynamical reasons. We discuss these next. 

Note that the T = vector breaks into two irreps 
for T > 0. This group theory need not contradict our 
intuition that at sufficiently low, but finite, temperature, 
the screening spectrum should be very similar to the T = 
hadron spectrum, and therefore the thermal and 
1 + should be nearly degenerate. When the screening 



1 The remaining symmetries are flavour symmetries. For staggered 
quarks the analysis is complicated |l4j ] by the mixing of flavour 
and spacetime symmetries. 
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mass, (i, is sufficiently small (i. e., p/T 3> 1), then the 
effect of the boundaries must be exponentially small. The 
approximate symmetry Oh is broken to D\ via terms that 
are exponentially small, leading to near-degeneracy of the 
components of T\ (the splitting can only be observed with 
exponentially large statistics). After all, the situation is 
not different from the T = theory defined with various 
different boundary conditions — they are equivalent as 
long as mL 3> 1 where L is the size of the box in which 
the theory is defined. 

However, this resolution raises another problem — can 
the masses of the two 0~ states, one obtained from the 
reduction of the PS and the other from the V be equal? 
If the mass of the latter becomes equal to that of the 1 + 
at low temperature, do we predict that the PS and V 
states must be degenerate, contrary to previous knowl- 
edge? The answer must be negative, the resolution being 
that the coming from the V cannot have overlap with 
the eigenvector of the transfer matrix with the smallest 
eigenvalue. This argument predicts that the screening 
masses in the V and PS channel (and similarly in the S 
and AV channels) must be different, although the naive 
transfer matrix argument gives a different result. 

The correlator identities 



C s (z) = -Cps(z) C v (z) = -C AV (z) 



(2) 



can be proven for overlap quarks by neglecting the effects 
of instantons. This shows parity doubling. As can be 
seen from Table U there is no group theoretical reason 
behind parity doubling. If there were, then that would 
force chiral symmetry at all T. 

Gauge field configurations which give rise to chiral zero 
modes of the overlap Dirac operator (called instantons) 
have interesting consequences. Since these modes are lo- 
calized, one recovers translation invariance by averaging 
over many such configurations. At high temperatures it 
is observed that these configurations are rather improb- 
able. Hence, to take their effects properly into account 
one must either use truly astronomical statistics or con- 
sider the zero momentum correlators, xr, i< e., correlators 
such as those in eq. ^ summed over 1 < z < N z . These 
have been analyzed in (2f| and various identities were 
checked. One should note that any configuration with a 
chiral zero mode breaks the identity Cs{z) = —Cps{z). 
Configurations with zero modes of both positive and neg- 
ative chirality break the identity Cv(z) = Cav( z )- Thus, 
generic ensembles of gauge configurations including in- 
stantons would realize only the irreps in Table [U Parity 
doubling, represented in eq. ([2]) involves the disappear- 
ance of chiral zero modes. It has been noted [25( that 
these do not disappear at the QCD phase transition but 
only gradually with increasing T. Hence parity doubling 
in the high temperature phase is an approximate state- 
ment, and extremely high statistics studies will be able 
to see the breaking of this approximate symmetry. 



III. FREE OVERLAP QUARKS 

The overlap Dirac operator (D) can be defined [24[ in 
terms of the Wilson-Dirac operator (D w ) by the relation 



d = i-d w (dId v 



-1/2 



(3) 



When the gauge fields are translation invariant, D w can 
be diagonalized in the Fourier basis. In free field theory, 
one can write 



Dw = c + i-ff+b^, where 

c = l + s-^sin 2 (p M /2), b p 



i Pp. ( 4 ) 



where s is the (negative) mass parameter, and, at finite T, 
the momenta are po — 2ir(n + 1 /2) /N t and pi = 2nn/N s . 
Since D'D = c 2 + b-b, and therefore diagonal, the overlap 
Dirac operator is 



D = a + i7 M /3 M , 
a = 1 



where 



(5) 



where \b\ 2 = b ■ b. Clearly the eigenvalues of D lie on the 
unit circle centered on the real axis at (1, 0). 

Writing the eigenvalues of D w as X(D W ) — rexp(i9), 
with r 2 = c 2 + \b\ 2 and 8 — tan _1 (|6|/c), it is easy to see 
that the eigenvalues of D are given by X(D) = pexp(i<j>) 
where p 2 = 2(1- cos 0) = 2(l-c/r) and = (-k-6)/2. It 
is easy to check that in the limit N t — > oo the minimum 
eigenvalue is obtained by taking po = iv/Nt and p% = 0. 
This gives 



7T 



P = 



(1 + s)N t 

- + o\ — 

2 \N t 



O 



(6) 
(7) 



The factor of 1 + s in the denominator is usually absorbed 
into a redefinition of the fermion fields, thus giving the 
correct Matsubara frequency in the continuum limit. 

A massive overlap operator, D(m), and the corre- 
sponding quark propagator, G(m), are defined by 



D(m) 
G{m) 
K{m) 



= m+ (1 - m/2)D, 



and 



= K-\m), 



where 



= [\-Dl2\D-\m), 



(8) 



and m is the bare quark mass in lattice units [241 ]. The 
massive overlap propagator can, therefore, be written in 
the form 



G(p; to) 
f(jp;m) 

h M (p;m) 



= f(p;m) +27^(^5 to), 



where 
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FIG. 1: Correlators computed in a free field theory of overlap 
quarks on lattice sizes 4 x 10 2 x 16, 6 x 14 2 x 24 8 x 20 2 x 32, 
12 x 30 2 x 48 and 16 x 40 2 x 64. The Af correlator coming 
from S (first panel), the E + correlator coming from V (second 
panel) and the correlator coming from V (third panel) 
are shown. The value of N T increases from top to bottom at 
zT = 2. The missing segments of the correlators in the last 
case correspond to the negative sign. 

As m — > 0, we recover the usual massless propagator. 

The screening correlator at external momentum q is 
then given by 

C r {q) = trG + 9; m)rG t (p; m)T. (10) 
p 

After performing the traces, one finds 

Cr(q)=±J2[f(p + q;m)f(p; m) + T^h^p + q)K(p)] , 
p 

(11) 

where the tensor r M „ is g^ for S and 2g^xg v \-gx\g^ for 
the A polarization of V. The tensors r M „ for the PS and 
AV satisfy the correlator identities for overlap quarks. 
Define CV to be the sum over three polarizations of the 



vector, then with the Euclidean metric one finds T^ v — 
-g^u - 2g zz . As a result, 

Cv(q) + C s (q) = -12^2 9z(P + l)9z{p), (12) 
p 

which is generically non-zero. 

In Figure [1] we show some of the results of a numer- 
ical computation of some of the meson-like correlators 
in free field theory on lattices of various sizes. One 
sees a gradual convergence of the results when plot- 
ted in terms of the scaling variables zT = az/N T and 
C(zT)/T 3 = a 3 C(az/N T )/N?. Several features of the 
free field theory are interesting enough that we list them 
explicitly. 

1. The correlator identities of eq. ^ are satisfied. 
This is expected, since free field theory has no topo- 
logical features which lead to localized zero modes. 
As noted before, these identities are stronger than 
the general feature expected from the group theo- 
retical analysis, since they allow us to equate the 

coming from PS to the ^4+ 1~ coming from the 
S, etc.. 

2. The correlators in the E + and irreps coming 
from the V have completely different behaviour. 
Consistent with this, the E~ and the Af corre- 
lators coming from the AV are quite different from 
each other. However, the E ± and the Af are iden- 
tical. 

3. The A[ correlator from the V changes sign at z ~ 
2/T. The point at which the change of sign occurs 
is almost independent of N z and N x , y . This is clear 
evidence that a three dimensional effective theory 
which describes these screening correlators cannot 
have reflection positivity. 

4. The Af coming from the V and AV are not degen- 
erate with the Af coming from the PS or the S, in 
agreement with the analysis of the previous section. 

5. The Af correlators from the PS and S are identical 
to the correlators from the V and AV. There is 
no group theoretical reason for this, and we shall 
examine later whether this holds in the interacting 
theory. 

IV. COMPUTATIONAL DETAILS 

Our numerical work was done in quenched QCD with 
overlap valence quarks. We have computed all 8 compo- 
nents of the correlators (one each for S and PS and three 
polarizations each for the V and AV). Extensive finite size 
scaling studies were reported earlier [25[ on lattices with 
a = 1/4T. We have extended this finite size scaling study 
here, but the main emphasis is on studying the variation 
with lattice spacing and taking the continuum limit. To 
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this end we have used a = 1/6T and 1/8T, at the tem- 
perature T = 2T C , where T c is the critical temperature 
for pure SU(3) gauge theory. 
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Lattice Size 


N 


6.0625 


4 x 10 2 x 16 


19 


6.3384 


6 x 14 2 x 24 


20 


6.55 


8 x 18 2 x 32 


26 


6.0625 


4 x 12 3 


23 


6.3384 


6 x 14 3 


16 



TABLE II: Coupling f3 for temperature 2T C , the lattice size, 
and the number of configurations (separated by 1000 sweeps) 
used for analysis (N). 

The computation of propagators using the overlap 
Dirac operator of eq. © needs a nested set of two matrix 
inversions for its evaluation (each step in the numerical 
inversion of D involves the inversion of D^Dw). This 
squaring of effort makes a study of QCD with dynamical 
overlap quarks very expensive. As before, we therefore 
chose to work with quenched overlap quarks. 

We generated quenched QCD configurations at tem- 
peratures of T/T c = 2 on 4 x 10 2 x 16, 6 x 14 2 x 24 
and 8 x 18 2 x 32 lattices (see Table ITT]) . The correspond- 
ing couplings are respectively f3 — 6.0625, 6.3384 and 
6.55. Note that these are the known critical couplings 
on N t = 8, 12 and 16 lattices. For the latter, no infinite 
volume extrapolation was available, unlike the first two. 
Our choice was motivated by the 2-loop /3-function, and 
consistency with the the finite volume results. In each 
case, the configurations were separated by 1000 sweeps 
of a Cabibbo-Marinari update. 

For the matrix M = D^D^,, and a given source vec- 
tor b, we computed y = M _1 ' 2 6 by using the Zolotarev 
algorithm [27| : 




where the coefficients q and di are computed with Jacobi 
elliptic functions once one has the order of approximation 
N and the ratio k = /Lt ma x/Mmin where /i max and ^ min 
are the boundaries of the domain where we apply the 
approximation. In our implementation of the algorithm 
(281 ] . we first compute the lowest and highest eigenvalues 
of M and then choose /x max and )i m i n so that the domain 
of applicability of the Zolotarev approximation is 10% 
larger than the domain spanned by the eigenvalues. The 
order No is defined by requiring a precision e/2 for the 
approximation of 1/ y/z in the entire domain. With these 
parameters, one calculates the approximation in eq. Q13p 
by a multishift CG-inversion at the precision e/2. 

On each gauge configuration wc computed G on 12 
point sources (3 colors and 4 spins) for 8 quark masses 
from m/T c =0.008 to 0.8 using a multi-mass inversion of 
D'D, The (negative) Wilson mass term in D w , i. e., l + s 



in eq. ([J}, which is an irrelevant regulator, was set to 1.8. 
The tolerance ranged from e = 10" 5 to 10" 7 in the inner 
CG and 10~ 3 to 10~ 5 in the outer CG, with most of the 
work done using the larger values of e. We remark on the 
choice in a later section. Typical No needed was about 
7-8. 

For most configurations we found that the spectrum 
of D* D starts well away from zero. However, very oc- 
casionally (one each for Nt = 4 and 8) we found a zero 
mode, i. e., an unpaired eigenvalue < 10 -2 . As discussed 
in the literature, such zero modes have to be subtracted 
off to obtain averages of physical quantities such as the 
correlation functions. We showed elsewhere that our de- 
termination of the eigenvectors and eigenvalues is precise 
enough that we can subtract out the zero mode contribu- 
tion in the chiral condensate. The resultant subtracted 
condensate (or correlators) is well within the statistical 
distribution of the measurement obtained in the sample 
without zero modes, at all the couplings and lattice sizes 
studied. We therefore chose not to include the zero mode 
contributions below. Furthermore, we show below results 
for our lowest quark mass m/T c =0.008, as the main fea- 
tures are essentially quark mass independent. 



V. RESULTS 
A. Symmetry realization 




FIG. 2: The A~[ and the E + correlators coming from the V 
correlator on an 8 x 18 2 x 32 lattice at T — 2T C . Also shown 
are the results of a free field theory (FFT) computation. 

We have checked that the sum of the S and PS cor- 
relators vanishes within errors in all our measurements, 
in agreement with the correlator identities for overlap 
quarks. For the V and AV correlators we checked that 
the difference of the x and y polarization correlators van- 
ish within errors, thus verifying that they lie in the same 
irrep of the symmetry group, as indicated by Table [IJ 
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FIG. 3: Screening correlators in the S and V channels com- 
pared with the corresponding ideal gas results on a 6 x 14 2 x 24 
lattice (first panel) and an 8 x 18 2 x 32 lattice (second panel) 
at T = 2T C . The lines are single cosh fits. 



In Figure [2] we show the screening correlator obtained 
from various polarizations of the V. The f-direction po- 
larization gives the and the x and y directions give 
the E + . In the figure we show the correlator from the 
sum of the x and y polarizations. We also show the cor- 
responding correlator obtained from a free field theory 
(FFT) computation. Note the close agreement between 
the E + correlator in the the quenched theory and the 
FFT results, over six orders of magnitude change in the 
correlator. We will examine this in more detail later. 

The disagreement between the FFT and the full theory 
for the correlator shown in Figure [5] requires further 
analysis. When the tolerance in the fermion inversion was 
decreased, the E + correlator changed only marginally, 
whereas the correlator changed sign at zT > 3/2 
on some gauge configurations. However, further work 
is required to check whether these correlators also turn 
negative as the tolerance is decreased further. 

In the past the sum of the three polarizations of the 
T = vector have been analyzed together. Strictly 
speaking, this is incorrect, since it mixes the E and 
irreps. However, the screening correlator in the latter ir- 
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FIG. 4: The ratio of the measured and ideal gas screening 
correlators in the S and V channels. 



rep is almost an order of magnitude smaller than that in 
the former. As a result, the screening masses extracted 
by the older method is a fairly accurate measure of the 
screening mass in the E channel. In fact, we find that 
the screening mass from the E + correlator at N t = 8 is 
afiE+ — 0.798(4) whereas that from the mixed represen- 
tation is afiv — 0.816(4). The two screening masses are 
compatible with each other at the 3ct level. 

With this caveat in mind, in the remainder of this 
paper we shall report results obtained with the older 
method of analysis, where the V correlator sums over 
all the polarizations. It has the advantage that the re- 
sults are directly comparable with earlier results, and is 
not far removed from the quantitatively correct results 
which would be obtained by following the theoretically 
correct method. 



B. Correlations as functions of spatial separation 

In Figure [3] we compare the measured correlation func- 
tions in the S and V channels. It is noteworthy that the 
V correlator agrees rather well with the ideal gas corre- 
lation functions from z = all the way to N z /2, when 
the correlator itself falls by 6-7 orders of magnitude. To 
make this spectacular point more forcefully, we have plot- 
ted in Figure [4] the ratio of the measured correlator with 
the ideal gas. It is clear that the V (E + ) correlator differs 
from the ideal gas by no more than 10-20% everywhere 
except (possibly) near the center of the lattice. 

The situation is somewhat different for the S (^4^) cor- 
relator. This agrees well with the ideal gas up to a dis- 
tance z ~ 1/T. At distances larger than this, the S chan- 
nel correlator falls significantly slower than the ideal gas 
result (Figure [3]). This is shown again in Figure [H where 
it becomes clear that the ratio of the correlators changes 
by a factor of about 5 when the correlators themselves 
have changed by about 6 orders of magnitude. 
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FIG. 5: Momentum space correlators at two different lattice 
spacings compared with the ideal gas results, normalized to 
fit the high-momentum data, for the S (first panel) and V 
(second panel) channels. 



C. Correlations as functions of momentum 

In some applications the momentum space correlator 
is a more direct input, giving emphasis to the small z 
region. Our lattice data can be exhibited in momentum 
space by a Fourier transformation of the separation z. 
The periodic boundary conditions give rise to a Brillouin 
zone structure as usual, and only momenta in the range 
< k z < wTNf are independent. 

We show the momentum space correlators in Figure 
[5j It turns out that the ideal gas result needs to be 
multiplied by a constant (independent of k z ) in order to 
describe the data well. However, after this normalization, 
which we choose to fix the correlator at k z = irTN t , the 
shapes of the correlators are in rough agreement with 
the ideal gas results, except at small k z . From Figure 
[5] another interesting feature which is visible is that the 
correlators at two different lattice spacings, a = 1/6T 
and 1/8T, after the scaling by a 3 , show no significant 
dependence on lattice spacing at k z = it /4a, -k /2a, 37r/4a 
and n/a. However, there are small differences visible at 
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FIG. 6: The ratio of the measured and ideal gas screening 
correlators in the S and V channels, on the 8.18 2 .32 lattice. 



low momenta, k z < 7rT, for both the S and V correlators. 
With decreasing a, these residual a-dependent effects go 
in the direction of slightly improving the agreement with 
free field theory for V, and of making matters slightly 
worse for S. 

The small disagreements with free field theory can be 
better exhibited by displaying the ratio of the correlator 
and its ideal gas counterpart- Cr(k z ) /Cf ree (k z ) . These 
ratios are shown in the V and S channels for the 8.18 2 .32 
lattice in Figured The ratio shows clear structure. Near 
the center of the Brillouin zone, k z w 47rT, the ratios go 
to a momentum independent constant. However, this 
constant is different from unity, and also different in the 
two channels. For k z < 2irT the ratio is clearly momen- 
tum dependent. In the S channel, the ratio decreases 
with k z whereas in the V channel it increases. To the 
best of our knowledge there is no prediction or under- 
standing of such structures in weak coupling theory. 



D. Screening masses 

One can extract screening masses from the correlators 
in two different ways. One is by fitting a single cosh to 
the long distance part of the screening correlator. The 
other method is to use the ratio C(z)/C(z + 1) and the 
assumption that at distance between z and z + 1 is de- 
scribed by a single cosh to extract a distance dependence 
effective mass, m(z). If there is a clear plateau in the 
effective mass as a function of z, which agrees with the 
fitted value, then one could reliably talk of a screening 
mass and its value. 

In Figure [7] we exhibit the result of such a test. It 
is clear that there is a stable plateau over almost half 
the available distances, and the local mass in this region 
(z > 1/T) agrees well with the results of a fit over the 
same range. This is quite different from earlier studies 
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FIG. 7: Effective masses extracted on the 8 x 1 8 2 x 32 lattice, 
compared with the results of the fits to a single cosh indicated 
by the horizontal lines. 



of staggered and Wilson quarks, where a stable plateau 
in the local masses was not observed. Note also the need 
for aspect ratios £ = N z /N s > 1 in order to stabilize the 
local masses at their asymptotic value. 

Notice the non-monotonicity of the local masses at the 
smallest distance. If the correlation function could be 
expressed in the spectral form — 

C(z) = M [exp(-miZ) + cxp(-rn l {iV 2 - z})] , (14) 
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with all the Ai > 0, i. e., the correlation function came 
from a transfer matrix satisfying reflection positivity, 
then the local masses would be a non-increasing func- 
tion of distance for z < N z /2. Clearly, therefore, the data 
shows that this condition is violated. Such behaviour has 
been observed with both staggered and Wilson quarks 
earlier. Note also the lack of convexity of the screening 
correlators in Figure[3l which is also a consequence of the 
violation of reflection positivity. 

In Figure [5] we plot the ratio of screening masses ex- 
tracted in quenched QCD to that in a free field theory 
of overlap quarks on the same lattice, [i/fi[ ree , as a func- 
tion of l/N? (x a 2 at fixed T — 2T C . In the V channel the 
ratio is consistent with unity. In the S channel, however, 
the ratio differs from unity at the 95% confidence level 
but is almost independent of the lattice spacing. 

We have shown results for two values of the aspect 
ratio, C — 1 an d 2. At the smallest lattice spacing we 
have only f = 2. It would be useful to extend these 
computations to larger £. If reflection positivity of the 
screening correlators were to hold, then one would expect 
the screening masses not to increase with £. However, in 
Figure [5] we see that this exactly what happens, again 
indicating that reflection positivity is violated. 

In view of our results, it is interesting to ask whether it 
is possible that at some larger value of N t there is a cross 
over to a regime where us is closer to /if ree . If such were 



FIG. 8: The ratio of the screening mass in quenched QCD 
to that in a theory of free overlap quarks on the same lattice 
for the S and V channels at T — 2T C . Results for two aspect 
ratios, f = 1 and 2, are shown. 



the case, then one should be able to identify the physics 
of such a cross over. There is only one such piece of 
physics, and that is the scale of non-locality of the overlap 
Dirac operator. This has been studied in [261 ] . where it 
was shown that the scale of non-locality varies smoothly 
towards zero as the lattice cutoff decreases. This allows 
us to rule out a cross over from the observed behaviour 
of us to trivial physics. 

A direct comparison of the screening masses in the S 
and V channels is shown in Figure [H From this it is 
clear that these two masses differ from each other at the 
95% confidence limit at the smallest lattice spacing that 
we use, and show no tendency to move towards equality. 
Since the ratio fis/fJ-v is within 5% of unity, an expla- 
nation within a weak coupling expansion seems possible. 
However, there is no such explanation available at this 
time. 
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FIG. 9: The ratio of the screening mass in the S and V chan- 
nels in quenched QCD are shown as a function of l/iV t 2 oc a 2 
for T — 2T C . Results for two aspect ratios, ( = 1 and 2, are 
shown. 

VI. SUMMARY 

We have performed a group theoretical analysis of the 
screening transfer matrix at finite temperature. Quark 
bilinear operators, in the C = +1 sector, which are clas- 
sified as S, PS, V and AV at zero temperature can be 
decomposed at finite temperature as shown in Table Q] 
The thermal state Aj~ is obtained in the decomposition of 
both the PS and V. Parity doubling is not a consequence 
of the group theory, but arises dynamically. Some iden- 
tities which can be proven for the topologically trivial 
sector of overlap quarks, eq. @, show parity doubling. 
We presented an argument that if the PS and V are not to 
be degenerate at low temperatures, then they cannot be 
degenerate at any temperature. The theory of free over- 
lap quarks is consistent with the converse, since it gives a 
degeneracy between the PS and V at all temperature. As 
a result, this theory cannot be an accurate representation 
of all aspects of QCD at finite temperature. 

We have demonstrated — 



1. The correlation function in the E + sector (V) is 
very close to that in free field theory. In the A+ 
sector (S) there is a significant difference at large 
distances, i. e., for z > 1/T. However, the magni- 
tude of such differences are small compared to the 
observed fall of the correlation function by 6 orders 
of magnitude over the distances involved. 

2. Both the S and V correlators are very close to the 
prediction of free field theory at distance less than 
T. 

3. The screening mass of the V is consistent with the 
free field theory of overlap quarks, whereas that of 
the S is lighter by about 5% at the 95% confidence 
limit. This behaviour persists into the continuum 
limit. 

4. The S and V correlators in momentum space agree 
with free field theory of overlap quarks, after an 
overall normalization, only near the center of the 
Brillouin zone. The part with k z < ttT differs from 
free field theory for both the S and the V. 

5. There is evidence that the screening transfer ma- 
trix does not have reflection positivity. This implies 
that screening in finite temperature QCD cannot 
be equivalent to the zero temperature physics of a 
temperature dependent Hamiltonian. This state- 
ment is not controversial, since screening phenom- 
ena at finite temperature crucially involve a mixed 
state density matrix, whereas correlations at zero 
temperature involve pure quantum states. 

None of these features depend upon topological non- 
triviality of the gauge field ensemble, since the one that 
has been used in these computations have vanishing 
topology. 

This work was funded by the Indo-French Centre for 
the Promotion of Advanced Research under its project 
number 3104-3. These computations were performed on 
the Cray XI of the Indian Lattice Gauge Theory Initia- 
tive (ILGTI) in TIFR, Mumbai. 
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